Theory of anisotropic s-wave superconductivity 
with point-node hke gap minima: 
analysis of (Y,Lu)Ni2B2C 
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Recent intensive experimental studies revealed that (Y,Lu)Ni2B2C is an anisotropic s-wave su- 
perconductor. In addition, its gap function possesses deep point minima, whose ratio of the gap 
anisotropy is more than 10. On the theoretical side, however, it is nontrivial to understand the 
origin of such a peculiar superconductivity. In the present paper, we propose a mechanism of 
the s-wave superconductivity with deep gap minima, based on the theoretical model where strong 
electron-phonon coupling as well as the moderate magnetic fluctuations coexist. By analyzing the 
strong coupling Eliashberg equation, we find that s-wave superconducting gap function owing to the 
electron-phonon coupling becomes highly anisotropic as the magnetic fluctuations increases. The 
set of model parameters for realizing the strong gap anisotropy in the present model will be ap- 
propriate for (Y,Lu)Ni2B2C. According to the present mechanism, (groups of) pair of gap minima 
appear at points on the Fermi surface which are connected by the nesting vector Q, in both cases of 
s-wave superconductors and non s-wave ones. We briefly discuss other superconductors with highly 
anisotropic gap function, e.g., PrOs4Sbi2 and A^ao.33Co02. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.40.-Fk 



I. INTRODUCTION 

A. superconductivity in RENi2B2C 

Since the boron-carbide superconductors RENi2B2C 
(RE=Lu, Y, Tm, Er, Ho, Dy) was discovered about a 
decade ago [1] , various experimental and theoretical stud- 
ies have been devoted to determine the symmetry or the 
mechanism of the superconductivity. The crystal struc- 
ture is body-centered tetragonal, consists of RE-C lay- 
ers separated by Ni2B2 sheets. Among them, LuNi2B2C 
and YNi2B2C show relatively high T^, 16. 5K and 15. 5K, 
respectively. They are non-magnetic metals till very 
low temperatures. On the other hand, in superconduct- 
ing Er and Ho compounds, /-electrons in Er and Ho 
show incommensurate magnetic long range orders with 
Qm ~ 27r(0.55/a, 0, 0), where a is the lattice spacing 
[2,3]. 

Various experimental results revealed that LuNi2B2C 
and YNi2B2C are s-wave superconductors, whose super- 
conducting (SC) gap functions are highly anisotropic. 
For example, the specific heat measurement in YNi2B2C 
below Tc tells that the SC gap in these compounds 
changes from a gap-less type to a full-gap one, by replac- 
ing Ni with Pt by 20% [4]. This result suggests that a 
s-wave superconductivity with deep point minima occurs 
in a pure compound, and its anisotropy is smeared out 
by impurities. Later, Izawa et al. measured the c-axis 
thermal conductivity of YNi2B2C in H rotated in various 
directions, and found that deep point minima in SC gap 
exist along [100] and [010] directions [5]. The estimated 
ratio of the anisotropy of the SC gap will be more than 
10. Recently, Watanabe et al. confirmed the point-node 



like SC gap function along [100] and [010] directions by 
the ultrasonic attenuation measurement [6]. 

On the theoretical side, however, the origin of such 
a highly anisotropic s-wave SC gap has not been un- 
derstood until now. It cannot be reproduced by solv- 
ing an Eliashberg equation even if one assume a highly 
anisotropic dispersion for conduction electrons and/or 
anisotropic electron-phonon (e-p) interactions. Thus, 
a reasonable theoretical model for the anisotropic s- 
wave SC gap with deep point minima, whose ratio of 
anisotropy is more than 10, is highly demanded for un- 
derstanding the superconductivity in (Y,Lu)Ni2B2C. 

In the present paper, we propose a mechanism of an 
anisotropic s-wave superconductivity with deep gap min- 
ima, in a system where strong e-p coupling and moder- 
ate antiferromagnetic (AF) fluctuations coexist. By solv- 
ing the strong coupling Eliashberg equation, we succeed 
in deriving a s-wave SC gap with deep point minima, 
whose ratio of anisotropy is more than 10, for a wider 
range of model parameters. The point- node like SC gap 
in (Y,Lu)Ni2B2C observed experimentally is reasonably 
explained by the present theory. The proposed mecha- 
nism of making deep gap minima due to the magnetic 
fluctuations is simple and general, so it will also be real- 
ized in various superconductors, including the unconven- 
tional superconductors. We study the case of the p-wave 
SC system in Appendix. 

By various experimental and theoretical studies for 
(Y,Lu)Ni2B2C, it is confirmed that both (i) strong e-p 
couplings and (ii) prominent A F fluctuations due to the 
nesting of the Fermi surface (FS) coexist in these com- 
pounds. Here, we explain the experimental and theoret- 
ical evidences of (i) and (ii) in more detail. 
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B. evidence for AF fluctuations 

NMR studies for (Y,Lu)Ni2B2C have been performed 
by several authors [7,8]. The spin-lattice relaxation ratio 
l/Ti of "B and that of ^^Y do not show Hebel-Slichter 
peaks. In the normal state, on the other hand, \/TiT 
of ^^B increases monotonously as temperature decreases, 
which suggests the enhancement of the AF fluctuations 
at lower temperatures. According to the spin fluctuation 
theory Hke the SCR theory [9], 1/TiT cx x%''^''^ , where 
d is the dimension of the system and XQ is the staggered 
susceptibility. In nearly AF metals, XQ shows the Curie 
Weiss temperature dependence. The observed 1/TiT of 
^^B can be fitted well by the above expression both for 
d=2 and 3. 

The strong AF fluctuations in (Y,Lu)Ni2B2C observed 
by NMR are expected to originate from the nesting of 
the FS: According to the LDA band calculations for 
LuNi2B2C [10 13], the FS possesses a nesting feature 
whose nesting vector is Q w 27r(0.5/a, 0), 27r(0, 0.5/a). 
The generalized susceptibility Xq(0) derived from the 
band structure given by the LDA study shows a peak 
at q ~ Q [10], which means that the RPA-type magnetic 
susceptibility, Xq(0) = Xq(0)/(1 ~ ^Xq(O))) lias a sharp 
maximum at q ~ Q. 

Another evidence of the nesting for q ~ Q is the mag- 
netic order of /-electrons in Er and Ho compomids, whose 
ordering vector is ~ 27r(0.55/a, 0, 0). This ordering 
is given by the RKKY interaction between /-electrons 
via the susceptibility of conduction electrons. Because 
all the band structures for RENi2B2C will be similar 
owing to the local nature of /-orbitals of RE, Xq(0) in 
(Y,Lu)Ni2B2C is expected to take the maximum value 
around q « Qm ~ Q- 



C. evidence for strong e-p coupling 

The density of states (DOS) at the Fermi energy in 
RENi2B2C obtained by the band calculation is 4.5 ^ 
4.8(eV cell)"-^ [11-13]. Because a cell contains two Ni 
atoms, the DOS per Ni is 2.25 ~ 2.4eV-^ This value of 
the DOS is much larger than that for high-Tc cuprates, 
which is about ^ 1.3(eV cell)~^ by the band calculation 
for La2Cu02, where a cell contains one Cu atom. 

Such a huge DOS in RENi2B2C suggests the strong 
e-p coupling. In addition, the light weight of B atom 
means the large frequency of the phonon mode, Wph. In 
fact, Wph ~ 300K is expected in these compounds [13]. 
The dimensionless e-p coupling constant A (or mass en- 
hancement factor due to e-p interaction) is defined as 
A = m*/mband — 1) where m* is the effective mass of an 
electron. Its value can be estimated by the temperature 
dependence of resistivity p using the Bloch-Gruneisen 
transport theory. The experimental value 0.4/Kf2cm/K 
in LuNi2B2C leads to A 2.6, although it might be too 



overestimated [11]. In addition, we stress that the promi- 
nent softening of the phonon dispersion at ii-ohs ~ Q is 
observed by neutron diffraction experiments [19]. This 
result strongly suggests the strong e-p coupling as well 
as the nesting of the FS with kobs ~ Q- Thus, the ob- 
served Kohn anomaly ensures the main character of the 
electronic properties in (Y,Lu)Ni2B2C. 

The estimations of the value of A using the thermody- 
namic measurements have been tried by many authors 
on the basis of the strong coupling Eliashberg equation. 
References [14] and [15] concluded that A = 1 ~ 1.2. 
On the other hand, A = 0.5 ^ 0.8 was deduced using 
the scanning tunnelling spectroscopy [16]. In addition, a 
moderate momentum dependence of (or FS dependence 
of ) A was inferred be several authors: Reference [17] ex- 
plained the experimental upper critical filed (-ffc2) by as- 
suming a two-band model with different A's (Amax = 0.8 
and Amin = 0.3). Later, more isotropic A was inferred by 
ref. [15]. More recently, Yamauchi et al. studied the mass 
enhancement factor A = ffldHvA/w-band — 1, where mband 
in the band mass given by the LDA study, and mdHvA is 
the cyclotron mass measured by the dHvA study, by as- 
suming that the mass enhancement is caused only by e-p 
couplings [18]. The obtained A's are 0.1~0.76 depending 
on the portion of FS's. 

However, previous works based on the strong coupling 
study could not reproduce the strongly anisotropic s- 
wave SC gap as observed in (Y,Lu)Ni2B2C, even if one 
assume a multi-band model with very different A's. In 
the present paper, we show the importance of the AF 
fluctuations to reproduce the deep point minima of the 
gap function found in (Y,Lu)Ni2B2C. 



D. s-wave or d-wave: theoretical viewpoint 

The existence of the prominent AF fluctuations in 
RENi2B2C suggest the possibility of the d-wave SC state, 
like in high-Tc superconductors. Fukazawa et al. per- 
formed the third order perturbation analysis with re- 
spect to [/ in a two-dimensional tight-binding Hubbard 
model for (Y,Lu)Ni2B2C, by neglecting the e-p couplings 
[28]. By solving the Eliashberg equation, they found that 
Tc ^ 15K for d-wave symmetry is realized by choosing 
a reasonable strength of U. However, the obtained 
might be overestimated because d-wave Tc is relatively 
low in the case of 3D systems in general. In fact, the 
FS of RENi2B2C given by band calculations possesses 
a three-dimensional (3D) structure, rather than a two- 
dimensional one. Experimentally, the anisotropy of the 
resistivity as well as that of Hc2 are small. In this sense, 
this compound is a three-dimensional superconductor. 

Next, we crudely estimate the s-wave due to the 
strong e-p coupling in RENi2B2C. By taking strong cou- 
pling effects into account, following BCS-McMillan [21] 
type expression for Tc valid for A ~ 0(1) would be ob- 
tained: 
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Tc = (a;ph/1.2) exp(-l/(A* - /x*)), 



(1) 



where A* = A/(l + A) and /i* is the Morel- Anderson 
pseudo-potential, which represents the reduction of 
owing to the pair-breaking effect by the Coulomb inter- 
action. By putting A = 1.5, Wph = 300K, and /i* = 0.15 
(a typical value) in eq.(l), we obtain Tc = 27K. If we put 
A = 1.0 instead, = 14K is obtained. 

In the following sections, we study the s-wave super- 
conductivity caused by the e-p coupling, under the influ- 
ence of prominent AF fluctuations. We propose an effec- 
tive model for (Y,Lu)Ni2B2C which reflect their char- 
acteristic properties well. By solving the strong cou- 
pling Eliashbcrg equation, we succeed in deriving the 
s-wave SC gap function with deep gap minima. The 
proposed mechanism of the strongly anisotropic SC gap, 
which is developed to explain the superconductivity in 
(Y,Lu)Ni2B2C, might be applicable to several unconven- 
tional superconductors found recently. 



II. THEORETICAL MODEL 
A. model with e-p coupling and AF fluctuations 

In this section, we explain the theoretical model used 
in the present work. In the present paper, we study the 
mechanism of the anisotropic s-wave superconductivity 
based on the simplified two-dimensional model under the 
influence of the magnetic fluctuations, which have not 
been studied previously. The proposed mechanism is ex- 
pected to occur in real compounds which have complex 
three-dimensional FS's. 

For the simplicity of the analysis, wc assume a two di- 
mensional isotropic Fermi surface as shown in Fig. 1(a). 
In addition, the density of states (DOS) at the Fermi 
level, iV(0), is assumed to be isotropic. Hereafter, we 
put = 2.5, which corresponds to n = 0.50 (quarter 
filling) if the area of the Brillouin zone is (27r)^ (square 
lattice). 

Next, we introduce the electron-electron interaction 

terms owing to the c-p interaction and the AF fluctu- 
ations. The latter originates from the (on-site) Coulomb 
interaction and the nesting of the FS. First, we repre- 
sent the interaction between electrons due to phonons, 
V'^^{ci>). In the case of the Einstein- type phonon. 
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where is the frequency of the Einstein phonon. g 
has a dimension of energy, and A = gN{0) gives the 
dimensionless coupling constant due to phonons. In 
(Y,Lu)Ni2B2C, Wph ~ 300K and gN{Q) ~ 1 is expected 
as mentioned in the previous section. 

We also introduce the interaction between electrons 
due to AF fluctuations, V^^{oj). In a Hubbard model 

with on-site interaction, U, it is given by ^§-Xq(w) 



within the RPA or fluctuation-exchange (FLEX) type 
approximation [22]. In the RPA or FLEX approxi- 
mation, the magnetic susceptibility Xq('^) is given by 
Xq('^) = Xq('^)/(1 - UXq{oj)), where Xq(t^) is the ir- 
reducible susceptibility. Here, we assume the following 
effective whose validity is assured in a lower energy 
region: 
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i + eirlQ-ql- 



iuj/ujsf 



model 1, 

model 2, 



(3) 

(4) 



where Q is the nesting vector where Xq=Q(0) takes the 
maximum value, ^af is the AF correlation length, and 
uigf represents the energy scale of the AF fluctuations. 
In both models, a has a dimension of energy, so aN{0) 
represents the dimensionless coupling constant for AF 
fluctuations. 

The model 1 is derived directly from Xq('*^) within 

the RPA or FLEX approximation, Xq(^) = Xq(^)/(1 ~ 
t^Xq(w)), by expanding Ux^ito) as C/x^(0) + ag^iplQ - 
qp - iasLo/oJsf + 0(|Q - ql^jW^), where as = 1 - 
Uxq{0) ^ 1 is the Stoner factor. Then, a is given by 

^3~Xq(0)- Thus, model 1 is has a correct functional 
form when |q — Q| ^ and uj ^ Wgf- The model 
1 had been studied intensively by Monthoux, Pines and 
their collaborators in the study of high-Tc superconduc- 
tors [23]. Within the SCR theory [9] or FLEX-type ap- 
proximation [22], cx 



and a 



^hich 



are observed experimentally in various high-Tc cuprates. 
Especially, Wgf > T (< T) in over-doped (under-doped) 
cuprates above the pseudo-gap temperatures, which is 
well reproduced by the FLEX approximation [24]. Ap- 
parently, u)st becomes small in the close vicinity of the 
magnetic quantum critical point (QCP). 

However, the model 1 might be unrealistic for |q — 
Ql » ^^p in that it has a lorenzian form with respect 
to |q — Ql, which decays too slow. Instead of introduc- 
ing a cutoff momentum Qc in the model 1, we introduce 
the model 2 which decays faster than lorenzian when 
|q— Ql ^ ^AF- "^^^ functional form of the model 2 will be 
justified if we assume the following expansion UXq{0) = 

t^x?j(o)+2asaFlQ-q|'+aieiFlQ-q|'+o(|Q-q|^), 

where we fixed the coefficient of the forth order term as 
•^S^AF- -'^s a matter of convenience, we put 2 in the nu- 
merator of cq.(4) to make the weights d^qV^^{0) for 
both models equivalent. In the present study, we treat a, 
^AF and Wgf as independent parameters both in model 1 
and in model 2. 

The momentum dependences of model 1 and 2 for 
a; = are shown in Fig. 2, where x = |q — Q|. x- 
depcndences of model 2 and model 3 (gaussian), both 
of which decay faster than lorenzian, looks similar. Al- 
though model 1 have been studied intensively in the 
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study of high-Tc superconductors [23], its supremacy is 
not apparent for the purpose of the present study. Note 
that a similar Tc for d-wave will be obtained by solving 
the Eliashberg equation in each models for the same pa- 
rameters in V^^{lo). In later sections, we study both 
model 1 and model 2, and analyze the latter model 
mainly. 

Finally, we point out that the isotropic FS used in the 
present study contradicts with the emergence of AF fluc- 
tuations, which should arise as a result of the nesting 
of the FS. This contradiction in the theoretical model 
will be discussed later. Although this simplification is 
unrealistic, it eliminates extrinsic complication in discus- 
sion and make the result clear. This simplification will 
not harm the generality of the the mechanism of the su- 
perconductivity with deep gap minima proposed in the 
present paper. 

B. origin of the deep SC gap minima 

Here, we explain qualitatively about the origin of the 

s-wave SC gap with deep minima. In the weak coupling 
BCS theory, where SC order parameter Ak is treated as 
energy-independent, it will be allowed to replace eqs.(2) 
and (3) with -ge{ujpi,-\uj\) and ¥^^^{0) -eiuj^i - \oj\), re- 
spectively [25]. As a result, Tc within the weak coupling 
BCS theory is given by the; highest temperature for the 
non-trivial solution Ak ^ of 

-E^k^-p(0)^tanh(^)^(a;.f-i;p), (5) 

where a singlet pairing is assumed. In the above equa- 
tion, Ep = ■^{cp — + A^, where ep is the dispersion 

of the conduction election and /i is the chemical poten- 
tial. The minus sign of the second term of eq.(5) can be 
interpreted as reflecting that {s ■ s') = —3/4 for a sin- 
glet Cooper pair. In the absence of the AF fluctuations, 
r° = 1.13a;phexp(— 1/A) for s-wave superconductivity is 
obtained by eq.(5) when A = gN{0) < 1 [25]. On the 
other hand, in the case of g = 0, the solution of eq.(5) 
is the dxy-'wavc like owing to the second term of eq.(5), 
^k^,ky = -A_fc^_fc„ = -Ak^^^ky = A_fc^,_fc^, as shown 
in Fig. 1(d). On the other hand, dx2_y2-type pairing is 
not favorable energetically when Q ~ (2A;f, 0), (0, 2fcF). 

Here we consider the s-wave solution of the gap func- 
tion Ak in eq.(5) owing to the strong e-p coupling. As the 
coupling constant for AF fluctuations, aN{0), increases, 
Ak with s-wave symmetry will has minima at (?k = n'K/2 
(n being a integer) by reflecting the negative contribution 
from the second term of eq.(5); see Fig. 1(b). Because of 
the simplicity of the mechanism, this theory of the s-wave 
SC gap with deep minima due to AF fluctuations is ex- 
pected to be general, independent of the fine shape of the 



FS. However, one has to check that the s-wave Tc, which 
will be smaller than T° (= Tc without AF fluctuations), 
is larger than the d-wave Tc caused by AF fluctuations. 

According to the present theory, each point minima is 
connected with others by Q. In the case of |Q| < 2kp, 
eight point minima will emerge (instead of four), except 
for IQI = \/2kF (see fig.l (c)). Except the number of 
point minima, the obtained results in the present study 
will qualitatively hold even if [Qj < 2fcF- Here, we note 
that the position of the point gap minima is equivalent 
to that of "hot spots" , where the quasiparticle damp- 
ing rate due to AF fluctuations take the (local) maxi- 
mum value. The concept of the hot spot is important 
to understand the anomalous transport phenomena in 
high-Tc cuprates [24]. When the AF fluctuations with 
Q — (Q,0),(0, Q) are strong enough, a d^y-typc super- 
conductivity as shown in fig.l (d) will be realized when 
Q < 2kp. As shown in fig. 1, the position of the nodes 
for the s-wave state do not coincide with that of the d^y- 
wave state, except for the case oi Q = 2kF. 

In the present paper, we will confirm that (i) a highly 
anisotropic s-wave SC gap observed in (Y,Lu)Ni2B2C 
(e.g., Amax/Amin > 10) cau be reafized with typical 
model parameters, and (ii) a relatively higher s-wave Tc 
(e.g., Tc > T"/2) can be realized under the condition that 
s-wave Tc is larger than the d-wave one. In the follow- 
ing sections, we solve the Eliashberg equation and de- 
rive the SC order parameter. The obtained region of 
the model parameters where (i) and (ii) are satisfied is 
wide enough and consistent with experimental situation 
in (Y,Lu)Ni2B2C. Hereafter, we study only the case of 
|Q] = 2fcF, which corresponds to Fig.l (b). 



III. STRONG COUPLING ELIASHBERG 
EQUATION 

In the preset section, we solve the Eliashberg equation 

for a singlet paring, and obtain the SC gap function and 
Tc for several sets of model parameters. Because the cou- 
pling constant A = gN{0) is of order 0(1) according to 
experiments, which corresponds to the strong coupling 
superconductor, we have to work on the strong coupling 
Eliashberg equation by taking the energy dependence of 
the gap function, Ai^{uj), into accoimt. 

The strong coupling Eliashberg equation for the 
present model within the one-loop approximation at zero 
temperature is given by [25,26] 



Ak(a;) = 




Re^M^=-X^_pico,z) (6) 

/■oo 1 

Ak-p(a;, z) = / dx-Im (V^'^ix) - V^Ip{x)) 
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ui + z + X — iS 



z ~ X + iS 



(7) 



for s-wave SC state, where gives the energy gap, 
which satisfies that = Ak(A^). To simpHfy the 
discussion, wc neglect the pair breaking effect by the 
Coulomb interaction U, that is, /i* = N{0)U/{1 + 
Af(0)[/ln(W^band/wph)) = is assumed [25,26]. (/x* is 
often referred as the Morel- Anderson pseudo-potential.) 
We also drop the impurity effect, that is, the clean limit 
case is studied. 

In solving the Eliashberg equation at T = 0, we will 
not use VP^{uj) given in eq.(2) because it gives an artifi- 
cial singularity of Ak(w) at lu = although it does not 
influence Tc as well as A^ badly. To escape the singular- 
ity, we use 



2 (a;-a;ph)2 + r2' 



(8) 



where F is a positive parameter which is much smaller 
than Wph- Equation (8) becomes gujpii(TT/2)6{uj — Wph) if 
r is infinitesimally small, which is equivalent to eq.(2). 
Hereafter, we put T = cjph/f 2. 

In eqs.(6) and (7), V|^p(x) is given in eq.(3) or eq.(4). 
Hereafter, we put |Q| = 2kF, which corresponds to Fig.l 
(b). Zk is the mass enhancement factor at w = 0, which 
is given hy = 1 — dT,i^{u))/du)\^^=o, where Sk(w) is 
the normal self-energy due to Vp^{u)) and V^Jp{w). It is 



expressed at T = as 

Z^ = l+gN{0) + N{0) 



f dnp 

Jfs 27r 



(9) 



Hereafter, we numerically solve the set of Eliashberg 

equations, eqs.(6), (7) and (9), imdcr the condition that 
the gap function has the (anisotropic) s-wave type sym- 
metry, i.e., A-k{u) has the four-fold rotational symmetry. 

In deriving eqs.(6)-(9), we performed the integration 
with respect to ek first by neglecting the imaginary part 
of the normal self-energy, which represents the quasipar- 
ticle damping rate. The validity of this approximation is 
apparently violated in high-Tc superconductors because 
the large ImSk(0)(::^ T) around the hot spots reduces 
A^(0). This effect decreases Tc much. The validity of this 
approximation for the present model will be discussed 
later. 

To obtain Tc both for s-wave and for d-wave, we solve 
the following Eliashberg equation at finite temperatures 
[26]: 

7riV(0) f dQp ApK„) 



~~ Jf 



Jfs 2n ^^2^ + AIM 

X (yP'^K - '^m) + V^IpicOn - COm)) , (10) 

(11) 



a-2«-i 



[i + aFlQ-q|']" + KIMf' 



(12) 



where w„ = 7rT(2n + 1), and a ~ 1(2) for model 
1(2). Tc is given by the highest temperature for the 
non-trivial solution of Ai^{ujn), under the constraint that 
^k:„ky = c ■ A_fc^,fe^ = c • Ak^,-ky = A_fc^ where 
c = 1 for the (extended) s-wave and c = — 1 for the Axy 
wave, respectively. In deriving eq.(lO), we used eq.(2), 
not eq.(8). 



IV. NUMERICAL SOLUTIONS 



A. comparison between model 1 and model 2 



Hereafter, we put the phonon parameters as A = 
gN{0) = 1.5 and Uph = 1-0; the latter corresponds to 
~ 300K experimentally. We stress that the main aim of 
this work is to present the new mechanism of the point- 
node like SC gap, not to reproduce the precise exper- 
imental value of Tc and thermodynamic measurements 
in (Y,Lu)Ni2B2C. Figure 3 shows A^^^x = 
Aj^ijj = A^^^Q, Tc for s-wave symmetry (s-Tc) and Tc for 
d-wave one (d-Tc) obtained in model 1 for = 400 and 
aiV(O) = 600. Wc see that s-wave SC state is real- 
ized when aiV(O) < 380. The condition A^^^/A^i„ > 10 
is realized when 300 < aN{0) < 380, although the re- 
alized Tc is less than one fifth of = 0.157, which is 
the transition temperature without AF fluctuations (i.e., 
aN{0) = 0). 

Next, we study the model 2: Figure 4 shows A^^^^, 
^min' s-Tc and d-Tc for w^f = 0.1 in model 2. We see 
that the condition A^g^^/A^j^ > 10 is realized for much 
wider range of aN{0) for = 50 ^ 200. In addition, 
the realized s-wave Tc is higher than the half of T^ . 

As a result, in both model 1 and model 2, a highly 
anisotropic SC gap is realized under the condition of 
s-Tc >d-Tc. However, the obtained s-wave Tc is much 
higher in model 2, and the condition for a strong 
anisotropy (say Aniax/Amin > 10) is much easier. The 
reason is that the slow tail of model 1 for q away from 
Q, which follows the lorenzian form, does not contribute 
to make SC gap minima, but work as a pair breaking 
effect like the momentum independent Coulomb interac- 
tion does, (see Fig. 2.) 

We note that the lorenzian type slow tail of model 1 
might not be natural because the peak of Xq(t^), which is 
created by the nesting of the FS, will decays faster once 
the nesting condition becomes ill as q is away from Q. In 
reality, Xq(0) for high-Tc cuprates obtained by the FLEX 
approximation, which is shown in Fig. 11 of ref. [24] for 
example, seems to decay faster than the lorenzian form. 
In that figure, we also see that Xq('^) takes an almost 
constant value (~ 0.5U~^) when |q — Q| ^ ^af^ 
will work as a pair breaking effect. As a result, the total 
Morel- Anderson pseudo-potential will be fXj^p+i^* 2/x*, 
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where fi* = N{0)U/{1 + A^(0)t/ln(lFband/^^ph). Here- 
after, we study the condition of reahzing the deep SC 
gap minima and analyze its property based on the model 
2 in more detail, by putting /x* = 0. The reader have 
to remind that the obtained Tc in later sections is over- 
estimated in that the pair breaking effect is neglected. 



B. gap functions in model 2 

Figure 5 shows Amax, Amin, s-Tc and d-Tc obtained in 
model 2 for ^|^p = 200. The o^sf-dependence of Amax.min 
and s-Tc is very weak, while that of d-Tc is strong. For 
Usf = 0.1, Amax/Amin > 10 is rcaHzcd for aN{0) > 280, 
and the condition s-Tc >d-Tc = 0.075 is satisfied when 
aN{0) < 490. 

Here, we discuss about the adequacy of parameters 
used in the present calculation. The relatively large 
value of ^AF used in the present analysis is consistent 
with the prominent AF fluctuations in (Y,Lu)Ni2B2C 
at lower temperatures [7,8]. Then, the smaller value 
of Usi is also expected because Wgf oc ^^p. Note 
that LOsf ~ T and oc are realized in an op- 
timum high-Tc cuprates, which is reproduced by the 
FLEX approximation [24]. Moreover, the coupling con- 
stant for AF fluctuations, aA'^(O), which corresponds to 
iV(0)(3t/V2)xQ(0) within the RPA or FLEX approxi- 
mation, will be smaller than that of high-Tc cuprates 
because U/Whand in (Y,Lu)Ni2B2C is smaller. As one 
can sec in Fig. 5, similar strongly anisotropic gaps are 
obtained for wider range of parameters when ^^p = 50 ~ 
400. 

Figure 6 shows the k-dcpcndcncc of Aj^ on the Fermi 
surface, for aN{0) = - 400. (Note that = 
ta,B~^ {kx / ky) .) Thus, the s-wave SC gap with deep point 
minima is realized for aN{0) > 300. It is similar to (sH-g)- 
wave SC gap function proposed in ref. [27], where the im- 
purity effect on the shape of A^. is analyzed. We stress 
that a model with almost the same strength of attractive 
interactions both for the s-wave and for the g-wave chan- 
nels, which was just introduced as an assumption in ref. 
[27], is derived microscopically in the present study. 

Figure 7 shows the e-dependence of Ak(e). Note that 
Ak(A^) = A^ gives the energy-gap of the SC order pa- 
rameter. We comment that F = Wph/6 in eq.(8) is as- 
sumed only in deriving Fig. 7, although F = a;ph/12 is 
used in other figures. (By this reason, A^ for aiV(O) = 
is smaller than 0.3 only in Fig. 7.) In Fig. 7, we note that 
Aniin(c) (for (?k = 0) increases as e increases from AjJ^j^^. 
Such an energy dependence of Amin(e) will stabilize the 
solution for the strongly anisotropic s-wave SC state in 
eq.(6). This effect, which is dropped within the weak 
coupling BCS theory, means the importance of the strong 
coupling effect for realizing the strongly anisotropic SC 
gap. 

At finite temperatures, the gap function Ak(e„) is ob- 
tained by solving the strong coupling Eliashberg equa- 



tion for Matsubara frequencies, eq (10). Then, the gap 
function for real e, Ak(e), is given by the numerical an- 
alytic continuation with high accuracy. Here, we study 
several physical quantities at finite temperatures using 
the obtained Ak(e): Figure 8 show the derived DOS for 
Caf = 200 and aN{0) = ~ 400, which is given by 



p(e)=^lmGk(u;-<5)A 



iV(0) / ?^Re 



Fs 2^ v^'-^ii^y 



(13) 



We put iV(0) = 1 in Fig. 8. It is shown that the finite 
DOS at lower energies emerges as the coupling constant 
for AF fluctuations, aA^(O), increases. We can also cal- 
culate the nuclear- lattice relaxation ratio (Ti) and the 
ultrasonic attenuation ratio (a). Figure 9 shows the ob- 
tained 1/TiT and a, which are given by 



R± = N{0) 



FS 



dOk 
"2^ 



dz ( 

az 



x2{[g^{z)r±[Mz)r), 
g]i{z) = J dekIniGk(-2 — iS)/iT = Re 

h{z) = j deklmFk(2: - i(5)/7r = Re 



Ak(.:) 
V^^-A^(^)' 



(14) 
(15) 

(16) 



where 1/TiT = i?+ and a = R-, respectively. We note 
that a in eq. (16) is derived by taking the average of the 
momentum of the absorbed phonon, which might be dif- 
ferent from the experimental situation. In Fig. 9, 1/TiT 
for aiV(O) = shows a large coherence peak below Tc. 
However, it is smaller than that given by the weak cou- 
pling BCS theory in the present model because the peak 
of the DOS is smeared out due to the strong coupling 
effect, that is, finite ImAk(e) at finite temperatures. The 
coherence peak is further suppressed as aN{0) increases. 

We comment that the coherence peak is further sup- 
pressed or vanishes if we take account of the effect of the 
finite quasiparticle damping rate, 7k = ImSk(— i(^), at 
finite temperatures, because finite 7k reduce the DOS at 
e = 0. However, this effect cannot be included in eq.(lO) 
because the following relation: 



&kGk(w„) = -TT , (17) 



TK^kF+AjK)' 



which is given by dropping the self-energy except for its 
tiin-lincar term, (— Zk -l- l)«w„, has been used in deriving 
eq.(lO). In cq.(17), Gk(i-^n) is the 2x2 Green function 
in Nambu representation, and is the Pauli matrix. To 
take the effect of 7k into account, which has been very im- 
portant in the case of high-Tc superconductors, we have 
to perform the ek-integration "numerically" instead of us- 
ing eq.(17). We further comment that the anisotropy of 
the dispersion for conduction electrons, which is totally 
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neglected in the present model, will decreases the coher- 
ence peak because a small k-dcpendcncc of Ak emerges 
even in the absence of the AF fluctuations. By these 
two reasons mentioned above, the tiny coherence peak 
for aN{0) > 300 may vanish after all. Note that a tiny 
coherence peak exists even in the case of d-wave within 
the BCS theory if the momentum dependence of A''(0) is 
totally neglected. 

Figure 9 also show the log(l/TiT)-logT plot for 
aiV(O) = ~ 400. This plot shows that the temperature 
dependence of 1/Ti at lower temperatures changes from 
an exponential behavior to a T'^-behavior for aN(0) > 
200, reflecting the finite DOS at lower energies around 
the hot spots. In a three-dimensional (3D) model, the 
relation 1/Ti oc corresponds to the line-node SC gap. 
If a point-node SC gap occurs in a 3D system owing to 
the present mechanism, 1/Ti oc will be observed. 

Before closing this section, we comment that the mass 
enhancement factor Zk is the function of the energy in 
the original Eliashberg equation. In the present work, 
however, we dropped the energy dependence for the sim- 
plicity of the analysis. It will be allowed for the purpose 
of the present study, although it is not for calculating the 
fine structure of the superconducting DOS. In Fig. 10, 
we show the solution of the Eliashberg equation given by 
eqs. (6) and (9) as well as the solution given by taking 
the energy-dependence of Zi^(lu) into account correctly. 
We see that the ratio of the gap anisotropy, A^^^^/A^jj^, 
is enlarged a bit, by performing a more accurate calcula- 
tion. In more detail, A^g^^ is enlarged about 20% whereas 
A^ij^ (at the hot spot) increases only shghtly, because 
the influence of the SC gap function on is little at the 
hot spot. We have also checked that the influence of the 
energy-dependence of Z^ on Tc is small. 



V. COMPARISON WITH THE WEAK 
COUPLING BCS THEORY 

In the previous section, we have shown that the s-wave 
superconductivity with deep gap minima is reproduced 
easily by solving the strong-coupling Eliashberg equa- 
tion for A = gN{0) = 1.5. That is, a large anisotropic 
ratio Amax/Amin > 10 is realized under the condition 
of s-Tc >d-Tc, by assuming reasonable model parame- 
ters. In this section, we study the same model within 
the weak-coupling BCS theory for a smaller value of A, 
and show that the strong-coupling analysis performed 
in previous sections is indispensable for reproducing the 
strongly anisotropic s-wave SC state. 

In the weak coupling BCS theory, energy dependence 
of the gap function as well as the normal self-energy are 
dropped. Thus, the corresponding Eliashberg equation is 
given in eq.(6) (or eq.(lO)) by (i) putting u; = (w„ = 0), 
(ii) replacing Ak(e) (Ak(en)) with A^, and (iii) putting 
= 1. The obtained Amin.max and s,d-Tc is shown 
in Fig.ll, which are denoted as "BCS". Here, we put 



A = gN{0) = 0.7 to justify the weak coupling treat- 
ment. We see that the s-T^ > dxy-T^ is realized only for 
aN{0) < 50, where the anisotropy of the s-wave SC gap 
is very weak. Thus, the strong anisotropic s-wave SC gap 
cannot be obtained by the BCS theory. 

Such a high dxy-Tc given by the "BCS" theory is 
suppressed drastically in the strong-coupling analysis, 
mainly due to the mass-enhancement factor; Zk > 1. 
To show this fact, wc also study the weak coupling BCS 
theory, by taking the "renormalization" by Zk into ac- 
count correctly. We call it the "R-BCS" theory, where the 
Eliashberg equation is given in cq.(6) (or eq.(lO)) by per- 
forming the simplification (i) and (ii) only. The obtained 
k-dependence of Z^ is shown in Fig. 12. We see that Z^ 
takes large values around the hot spots, so it will mainly 
reduce the effect of the AF fluctuations in the Eliashberg 
equation. The obtained result for gN{0) = 1.0 by the 
R-BCS theory is shown in Fig. 11. We 
is strongly reduced by the factor Zk, so s-Tc >dxy-Tc is 
realized till much larger vahie of aA^(O): aA^(O) < 240. 
Thus, Amax/Amin ^ 5 cau bc realized within the R-BCS 
theory. 

As a result, the renormalization effect by Zk is im- 
portant to realize the highly anisotropic s-wave SC gap. 
However, R-BCS theory is still insufficient for a quan- 
titative study, so one have to work on the strong cou- 
pling Eliashberg equation, eq. (6) or eq. (10), for 
gN{0) ~ 0(1). 

Finally, we comment that the DOS around the hot 
spots will decrease appreciably owing to the finite 7k be- 
cause 7k takes a large value at the hot spots when the 
AF fluctuations are strong, as .Zk does. This effect was 
found to reduce the Tc prominently in the case of high- 
Tc superconductors. To take this effect into account, we 
have to perform the Ck-integration "numerically" instead 
of using eq.(17). In this respect, the obtained dxy-Tc in 
the present work may be over-estimated. 



VI. DISCUSSIONS 

A. condition for the strongly anisotropic s-wave SC 
gap in 3D systems 

In the present paper, we proposed a mechanism of the 
s-wave superconductivity with deep gap minima in the 
presence of AF fluctuations. We found that a highly 
anisotropic SC gap observed in (Y,Lu)Ni2B2C can be re- 
produced with reasonable model parameters. Here, we 
summarize the obtained condition for realizing the deep 
SC gap minima: 

(i) To make A^g^x/^min ^ 1; a strong e-p coupling (e.g., 
A = gN{0) > 1) as wefi as larger (e.g., > 50) 
is required. The "radius of the minima" in the SC gap 
will be ~ ^^p. The strong coupling Eliashberg equation 
should be solved for a reliable analysis. 
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(ii) To make s-Tc >d-Tc; in addition to a strong c-p cou- 
pling, a smaller coupling constant between electrons and 
AF fluctuations (e.g., aiV(0)/^^F < 1) ^ ^^11 as a small 
energy scale of AF fluctuations (e.g., Wgf '^ph) is re- 
quired. In addition, the nesting area of the FS should be 
small. This condition will be easily satisfied in a three- 
dimensional system with several FS's. This fact will be 
discussed later in more detail. 

Here, the shape of the FS is chosen to be circle. The 
most essential effect of the anisotropy as well as the di- 
mensionality of the FS will appear in d-T^., because it is 
enhanced when the nesting of the FS which is consistent 
with the main AF fluctuations exists, while s-Tc and s- 
wave gap function will be rather insensitive. Thus, the 
d-Tc obtained here may have only a qualitative meaning. 
Nonetheless, the realized d-Tc will remain small if one 
analyze a three dimensional (3D) system. In fact, refs. 
[28-31] discuss that the d-Tc due to AF flTictuations in 3D 
systems remains low because the nesting area in 3D FS 
is small in general. In addition, the region of parameters 
for d-Tc is also very restricted in 3D systems. According 
to a band calculation for YNi2B2C, the shape of the FS 
is 3D like, and the nesting area of the FS is only 4.3%, 
which is supported experimentally by a tiny change of 
the resistivity at T = Tn [13]. 

As a result, d-Tc will not be high in the real three 
dimensional FS for (Y,Lu)Ni2B2C, so the mechanism of 
s-wave SC gap with deep minima proposed in the present 
paper will be realized in (Y,Lu)Ni2B2C without difS- 
culty. In general, the weight of the nesting area in the 
FS decreases further if several FS's exists like in heavy 
fermion systems. Although the present study is based 
on a simplified two dimensional model, the mechanism 
of the anisotropic s-wave SC gap due to AF fluctuations 
proposed in the present work will be universal. We also 
not that the effect the finite 7k — ImEk(— i<5) at T = Tc, 
which will reduce d-Tc as mentioned in the previous sec- 
tion, is not taken into account in the present analysis. In 
this sense, the obtained d-Tc is overestimated. 

Finally, we shortly discuss the location of the gap min- 
ima in 3D systems. Figure 13 shows the expected loca- 
tion of the gap minima for (a)|Q| = 2fcF and (b)|Q| < 
2kF (Q II X in both cases) in the case of the spherical 
FS whose radius is fcp. In case (a), point minima occur 
at two points on the .T-axis. In case (b), on the other 
hand, gap minima will make two small circles. There- 
fore, gap minima will be short lines or small circles in 
real 3D systems with anisotropic FS's. Obviously, the 
depth of the gap minima should depend on the position: 
The deepest point will be on the portion of the FS where 
A or the density of states at k, lmGk{—i5)/7T, takes rela- 
tively a smaller value. In such a case, the deepest point 
of the gap may cause "point-node like behaviors" in var- 
ious thermodynamical quantities at lower temperatures, 
even in case (b) (|Q| < 2A;f). 

For a more detailed study, a strong coupling study 
based on a microscopic Hamiltonian with 3D FS's is 
highly demanded. For example, an analysis of a Holstein- 



Hubbard model by the FLEX approximation will be 
fruitful. 



B. possibility of the (s-|-id)-wave superconductivity: 
breakdown of the time reversal symmetry 

In figs. 4 and 5, we obtained s-Tc and d-Tc by con- 
sidering the transition from the normal state. In the 
case of s-Tc >d-Tc (d-Tc >s-Tc), a pure s(d)-wave SC 
state will be realized whole the temperature region be- 
low s(d)-Tc because the established s(d)-wave SC gap will 
prevent the emergence of the d(s)-wavc gap. However, it 
is an interesting question whether s-wave SC gap and 
d-wave one coexist or not when d-Tc -^s-Tc below min{d- 
Tc, s-Tc}. In this case, the gap fimction will be given by 

= Ak3'^ -I- e*'^A^ , where wc choose the U{1) gauge 
to make both A^" (cx 1) and Aj^" (cx kxky) real. Then, 
gi^ = i is rcquirc;cl to make; the gain of the condensation 
energy largest, at least within the weak-coupling theory, 
because (A^)fs takes the maximum value then. The 
(,s + Z(i)-wave supcrconduc;tivity is a very interesting state 
in that the time reversal symmetry (TRS) is broken. In 
a later publication, we will study the nature of the super- 
conducting state whole the temperature region in more 
detafl [32]. 

Here, we notice that some of the deep minima on a 
s-wave SC gap caused by the mechanism proposed in the 
present work might be filled to some extent once the sys- 
tem changes from the s-wave SC state to the (s -|- id)- 
wave state as the temperature decreases. In the present 
model, in fact, the position of the "point minima" in 
the anisotropic s-wave state is different from that of the 
dx2-y2-wa.ve state except for the case of |Q| = 2A;f, as 
shown in Fig. 1. We expect that a pure s-wave state is 
realized in (Y,Lu)Ni2B2C because the nesting area of the 
FS is very small, which suggests that s — Tc <C d — Tc as 
discussed before. 



C. consideration on other mciterials: Sr2R-u045 
PrOs4Sbi2, Nao.33Co02 

First, we stress that the proposed mechanism here 
can make deep SC gap minima also in a unconventional 
(p-wave or d-wave) superconductors besides the original 
nodes: When a unconventional superconductivity is real- 
ized by a Kohn-Luttinger type mechanism, which is free 
from the concept of the QCP, some magnetic fiuctua- 
tions owing to the nesting of the FS will make "deep gap 
minima" on the SC gap function, if their energy scale is 
too low to produce a different type unconventional super- 
conductivity. Many heavy fermions and transition metal 
oxides superconductors have been reproduced properly in 
terms of the Kohn-Luttinger type mechanism, using the 
perturbation theory with respect to U [33] . In Appendix, 
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wc study a model which shows the p-wave superconduc- 
tivity with full gap, and find that the deep gap minima 
emerges on the SC gap function as the AF fluctuations 
increases. 

In a similar context, we would like to explain the theo- 
retical study for Sr2Ru04 by Nomura and Yamada [34]: 
Sr2Ru04 shows a spin-triplet SC state at Tc = 1.5K. 
Based on the third-order and forth-order perturbation 
theory, they found that the p-wavc superconductivity oc- 
curs in Sr2Ru04. According to their analysis, a large 
and isotropic (chiral p-wave) SC gap opens on the 7- 
FS, whereas small and strongly anisotropic gaps occur 
on a,/3-FS's. The obtained SC gap function is consis- 
tent with the specific heat measurement below T^. It is 
also consistent with the recent heat conductivity mea- 
surement under the rotatable magnetic field [35]. We 
note that the almost gap-less SC gap realized in Sr2Ru04 
originates from the cancellation of the pairing interac- 
tion. 

PrOs4Sbi2 is a Pr-filled Skuttcruditc superconducting 
compound, with = 1.85K [36,37]. In PrOs4Sbi2, con- 
siderable magnetic (dynamical) fluctuations are observed 
by /x-SR measurement [38] and by neutron diffraction 
measurement below 4K [39]. Moreover, e-p coupling will 
be also strong in PrOs4Sbi2 because various Pr-filled 
Skutterudite compounds are conventional s-wave super- 
conductors with relatively high T^s. 

The symmetry of the SC gap is not determined now. 
According to the heat conductivity measurement under 
the rotatable magnetic field, the SC gap has four point 
nodes along [100] and [010] axes in a higher magnetic-field 
phase, whereas two point nodes along [010] axis disappear 
in a lower magnetic-field phase [41]. One may expect 
that the strongly anisotropic s-wavc SC state is realized 
in PrOs4Sbi2, considering that the small energy scale 
of the AF fluctuations observed by neutron diffraction 
(wsf ~ 0.5meV) [39] is appropriate for making the deep 
SC gap minima as shown in the present paper. More- 
over, a recent /i-SR measurement suggests that the TRS 
is broken in the SC state of PrOs4Sbi2 [38]. It may be a 
(s-l-id)-wave SC state [40] , which could be reproduced by 
the mechanism proposed in the present work as discussed 
in the previous subsection. 

Nao.33Co02 is a triangular lattice cobalt oxide super- 
conductor, with Tc = 4.5K. The symmetry of the SC 
gap is also under debate now. According to NMR/NQR 
measurements, 1/T\T in a sample with magnetic fluctu- 
ations shows no coherence peak, and 1/Ti cx is ob- 
served below Tc [42,43]. On the other hand, 1/TiT in a 
sample with less magnetic fluctuations shows a tiny co- 
herence peak [44] . The reduction ratio of Tc due to impu- 
rities, —dTc/dx {x being the concentration of impurities), 
seems to be too small as a unconventional superconduc- 
tor. These measurements might be able to be explained 
as a s-wave superconductor with deep gap minima caused 
by the mechanism proposed in the present paper, as is 
the case with (Y,Lu)Ni2B2C. 

In both PrOs4Sbi2 and Nao.33Co02, more experimen- 



tal and theoretical studies are required to determine the 
symmetry and the mechanism of the superconductivity. 
In future, the present study will offer a hint to under- 
stand various SC compounds in the presence of both e-p 
couplings and magnetic fluctuations. 



D. summary 



In summary, we studied the influence of the AF fluc- 
tuations on the SC gap function Ak(e) with s-wave sym- 
metry owing to the strong e-p coupling, by analyzing the 
strong coupling Eliashberg equation. We confirmed that 
deep SC gap minima emerge in Ak(e) as the AF fluc- 
tuation increases. The condition for the realization of 
strong anisotropy (say A^jj^^/Ajjjjjj > 10) is studied in 
detail. We stress again that the main aim of this work 
is to present the new mechanism of the point-node like 
SC gap, not to reproduce the precise value of Tc and 
thermodynamic measurements in (Y,Lu)Ni2B2C. 

According to the present mechanism, (groups of) pair 
of gap minima appear at points on the Fermi surface 
which are connected by the nesting vector Q, in both 
cases of s-wave superconductors and non s-wave ones. 
Their position is equal to the hot spots, where the quasi- 
particle damping rate 7k caused by the magnetic fluc- 
tuations takes the (local) maximum value. The present 
mechanism is expected to reproduce the point-node like 
SC gap in (Y,Lu)Ni2B2C, as well as the direction of the 
point minima. It is a future problem to determine the 
precise positions of the gap minima on the FS's by taking 
account of the realistic shape of FS's. 

To realize a higher s-Tc (»d-Tc), it is desired that the 
position of the point minima (hot spots) is away from the 
van-Hove singular point, because A is expected to take a 
larger value there. In case there are several hot spots, the 
deepest point gap minima will appear on the FS where 
A takes the smallest value. Deep SC gap minima in the 
SC gap due to AF fluctuations proposed in the present 
work are expected to emerge even in unconventional su- 
perconductors, as explained in Appendix. In addition, 
recently found exotic superconductors, PrOs4Sbi2 and 
Nao.33Co02, were briefly discussed from the viewpoint 
of the present study. 

The author is grateful to Y. Matsuda, K. Izawa and 
J. P. Prison for useful discussions about experiments 
on (Y,Lu)Ni2B2C. He is also grateful to M. Sato, H. 
Yoshimura and K. Ishida for valuable comments on 
NMR/NQR. For useful discussions on theory, he is thank- 
ful to K. Yamada, T. Sato, H. Kohno, H. Ikeda and T. 
Nomura. 
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APPENDIX A: IN THE CASE OF THE P-WAVE 
SUPERCONDUCTOR 

As discussed in §VI, the deep SC gap minima at the 
hot spots owing to the AF fluctuations proposed in the 
present paper will also be realized in the case of p- or 
d-wave superconductor. To confirm this fact, wc study 
a p-wave superconducting systems with AF fluctuations, 
and see that the deep SC gap minima emerge as the AF 
fluctuations increases. In general, the odd-parity SC gap 
function in the Nambu representation is given by 



(Al) 



where dk is the d-vcctor. The quasiparticle energy gap 
is given by y/d^ ■ dk. 

Here, we introduce the following paring potential for 
the p-wave channel in two dimension: 



2cos(0k-ek'), 



(A2) 



whose diagrammatic expression is shown in Fig. 14. In 
the presence of the AF fluctuations, the system has the 
tetragonal symmetry [Dah)- In this case, six p-wave pair- 
ing states with different symmetries give the same quasi- 
particle gap without nodes, which is expected to optimize 
the condensation energy at least within the weak coupling 
theory [45]. In each pairing state, dk(w) is composed of 

(^bS^ {ui)^ A|^^''(a;)), which changes under the point group 

operations as {kx,ky). We promise that both A^^'(O) 

and A^^^(O) are real. For example, the chiral p-wave 

state is given by dk(cc)) = z{A'^\lj) ± iA'^\uj)). Then, 
the strong coupling Eliashberg equation for the p-wave 
symmetry with full gap is given by 



aL^'^K) 



7riV(0) 



Jfs 27r 



^p-ch 



-1, 



where J = 1,2, Dl{io) ^ {A^\u:)Y + {A'^\u)}\ 
and Vj^p(wra) represents the interaction due to AF 
fluctuations introduced in cq.(12). Note that the fac- 
tor — 1/3 in front of V-^^^ comes from the fact that 
(s • s*)g_2 = • s'^g^Q = 1/4. Because of the 

relation cos((?k ^ ^^-k) ~ —1, the expression for Zk in 
eq.(A3) is same as that in eq.(9). 

Here, we study eq.(A3) numerically, by putting g = 
Up = 1.0. The quasiparticle gap is given by the relation 
A^ = |£)k(A^)|. The upper panel in Fig. 15 shows the 
obtained A^^, = Ag^^^ and A^i„ = A°^^^^^ for the p- 
wave SC state at zero temperatures, as functions of the 
AF couphng constant aN{0). (Note that | Ak| is constant 
when aN{0) = 0.) The lower panel shows the obtained 



Tc's both for p-wave and d-wave. As shown in Fig. 15, 
deep minima emerges in the p-wave SC gap owing to 
the AF fluctuations. Quantitatively speaking, however, 
the obtained ratio Amax/Amim in the present model is 
smaller than that in Fig. 5 at a fixed aN{0), because of 
the factor 1/3 in front of Vaf in eq.(A3). 
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FIG. 1. (a) Fermi surface of the present study. We put 
fcp =, which corresponds to the quarter-filled band in a square 
lattice system, (b) anisotropic s-wave SC gap in the presence 
of the AF fluctuations; Q = [Ikp, 0), (0, 2feF). (c) anisotropic 
s-wave SC gap in the case of Q = (^fcF,0), (0, v/2fcF). 
(d) dxy-wave SC gap caused by the AF fluctuations; 
Q = (Q,0),(0,Q). 







model 1 


^=1 


model2 /'n 


''A 




models / ' 

,/ 




// 


V 

\\ 




\\ 






"/ 


/ 


'^.--t!! . 



-4-2 2 



FIG. 2. 14% (0) for model 1 and model 2, with 
o = 1 and ^AF ~ 1. Model 3 is defined as 
^x+q(0) = V^a • exp(-5ip|x|^). 
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FIG. 4. A^in and s,d-Tc obtained in model 2 for 
cjsf = 0.1. FIG. 6. 6k-dependence of AjJ. obtained in model 2. 
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FIG. 7. e-depcndcnce of Amin,max(£) obtained in model 2. 
Here we put F = uJaf/G- 
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FIG. 9. (a) Obtained nuclear-lattice relaxation ratio 1/T\T 
and ultrasonic attenuation ratio a. (b) 1/Ti is proportional 
to at lower temperatures for aN{0) > 200, reflecting the 
deep minima in SC gap. 1/Ti oc should be observed in a 
three-dimensional model. 




FIG. 8. e-dependence of p(e) obtained in model 2. We put 
Ar(0) = 1. 



FIG. 10. 6'k-dependence of obtained by solving Eliash- 
berg equation when aN{0) — 300, (i) by dropping the en- 
ergy-dependence of Zk (dashed line) and (ii) by taking the 
the energy-dependence of Zk{u)) into account correctly (full 
line). 
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FIG. 11. Obtained and s.d-Tc by BCS theory {Zk = 1) 
and by R-BCS theory {Zk > 1), respectively. 
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FIG. 12. 6'k-dependence of Zy, for A = 1.5. 
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FIG. 13. Dlustration of the location of the gap minima for 
(a)iQi = 2kF and (b)|Q| < 2kF (B || x), in the case of the 
spherical FS. 
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FIG. 14. Paring interaction for the p-wave channel. 
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FIG. 15. Obtained quasiparticle gap function for the 
full-gap p-wave SC state and Tc's for both p- and d-wave so- 
lutions, as functions of aN{0). 
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